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We investigate the deformation of wavefronts of sound waves in rotating Bose-Einstein conden¬ 
sates. In irrational fluid flows Berry et al. identified this kind of deformation as the hydrodynamic 
analogue of the Aharonov-Bohm effect. We study this effect in Bose-Einstein condensates and 
obtain the Aharonov-Bohm phase shift at all wavelengths. We show that this deformation of wave 
fronts is seen in both phase and density fluctuations. Eor wavelengths larger than the healing length, 
the phase fluctuations experience a phase shift of the order of 27r times the winding number. We 
also consider lattices of vortices. If the angular momentum of the vortices are aligned, the total 
phase shift is times the number of vortices in the condensate. Because of this behaviour the hy¬ 
drodynamic Aharonov-Bohm can be utilized as a probe for quantum vorticity, whose experimental 
realization could offers an alternative route to investigate quantum turbulence in the laboratory. 
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In 1980 Berry, Chambers, Large and Upstill [T] realized 
that surface waves scattered on a vortex pick up a geo¬ 
metric phase shift. This geometric phase shift turns out 
to be the exact analogue of the Aharonov-Bohm phase 
shift [2]. Moreover, the absence of gauge symmetry in the 
hydrodynamic analogue allows for new features that are 
unobservable in the quantum mechanical problem. In a 
fluid, the wave undergoes a wavefront dislocation [3] that 
is, entirely new wavefronts come out of the vortex core. 
This hydrodynamic Aharonov-Bohm effect has been con¬ 
sidered by several authors, for water waves Hi] or for 
superfiuid helium HID. The main difference between 
classical and quantum fluids, is that the vortex circula¬ 
tion is quantized. According to Berry’s classification [9], 
there are four types of Aharonov-Bohm effect, depending 
whether the wave and/or the flux (here the vortex circu¬ 
lation) is quantum or classical. Here we investigate the 
Aharonov-Bohm effect on sound waves (classical) scat¬ 
tered on Bose-Einstein condensates vortices (quantum). 

Wavefront dislocations have been used theoretically 
and experimentally in Bose-Einstein condensates to 
detect vortices m and their solitonic version El. This 
dislocation arises from a geometric phase in the conden¬ 
sate wave function. Here instead, we shall investigate 
how vortices affect plane wave perturbations. The aim 
of the present paper is twofold. We present a detailed 
analysis of the analogue Aharonov-Bohm effect, which 
includes dispersive effects, and compare it with the 
dislocation of the condensate wave function itself. The 
hydrodynamic Aharonov-Bohm effect in a Bose-Einstein 
condensate has so far not been considered. As we show, 
for wavelengths of the order of the healing length, 
the induced phase shift is of order unity and should 
therefore be experimentally measurable. This effect also 
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cumulates when considering a lattice of vortices, which 
is the stable configuration of rapidly rotating conden¬ 
sates E|. This allows us to probe vortices using sound 
waves, and obtain not only the total angular momentum, 
but also its orientation, and hence offers an alternative 
method to explore quantum turbulence in the laboratory. 

Quantum vorticity. — In the Bogoliubov approxima¬ 
tion, a Bose-Einstein condensate is described by its con¬ 
densate wave function T m- Its modulus gives the atom 
density, and its phase the velocity field, that is 

T(t, 7) — \/n(7^) exp {—ijit + i0( r^)). (I) 

Here /i is the chemical potential, and n{'r) the local num¬ 
ber density. In umts where h = 1, the velocity field 
is given hy If = V0/M, where M is the mass of the 
atoms of the condensate. In a condensate, the flow is al¬ 
ways irrotational, as it is by definition a gradient. How¬ 
ever, wherever the density vanishes, the phase becomes 
ill-defined, and the velocity field can acquire a nonzero 
circulation. This give rise to vortices in the condensate 
(see [12] for a review). If a single vortex is present, the 
condensate wave function reads 

r) = (2) 

where we used cylindrical coordinates {r,0,z). Because 
the condensate wave function must be single valued, w 
is an integer. This means the circulation of the velocity 
field is quantized and given by 

= 2ecre, (3) 

where c is the speed of sound in the condensate and 
^ = (2Mc)“^ the healing length. Around a vortex, the 
density is essentially constant n(r) ^ no, except within a 
region of size ^ which corresponds to the core of the 
vortex [T4| . Inside that region, the density decreases until 
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it vanishes at r = 0, the center of the vortex. Moreover, 
in a real condensate, quantum depletion further dresses 
the flow near the core of the vortex. The non-condensed 
part of the atoms induces a flow that is no longer re¬ 
stricted to be curl-free, and one can show that close to 
the center, the flow essentially becomes that of a solid ro¬ 
tation [E]. In the sequel, we shall consider perturbation 
far from the vortex core, so that we can safely approxi¬ 
mate the density by n(r) ^ no and the velocity flow by 
its superfluid part 


V = -60, 

r 


(4) 


which is the expression for an irrotational vortex (e^ is 
the unit angular vector). When two such condensates 
interfere, the pattern shows a wavefront dislocation on 
each vortex core [IQlIIl]. The number of new wave fronts 
is precisely given hy w. As we shall see, perturbations 
around the condensate undergo a similar effect. How¬ 
ever, depending on the wavelength, the geometric phase 
will smoothly interpolate between 0 , and 27rr6, i.e. w 
dislocated fronts. 


Sound waves .— It is well-known that excited states of 
a condensate can be obtained from the eigen-modes of 
the linearized Gross-Pitaevskii equation m- For this we 
assume a condensate wave function of the form 


^'0, r) = ?/)o0, r)(i + ut^(r)e + (5) 

where ^po is given by Eq. The fluctuations modes 
Vuj then obey the Bogoliubov-deGennes equation 

oouuj = Uuj + Mc^Vuj, ( 6 a) 

• V} + + Mc^u^. ( 6 b) 

In the absence of a flow, the modes can be decomposed 
in plane waves = Uke^^' ^ of momentum k. The 

above equations then impose the Bogoliubov dispersion 
relation 


ujl = c\e^ek^)=F{e). (7) 

We point out that we do not decompose u^j and v^j in 
circular waves of fixed azimutal number (eigen-modes of 
ide) as is usually done [161 [17]. This is because we want 
to see how plane wave excitations are affected by the 
vortex. We shall now proceed in two steps. We first 
solve the Eq. (§ to obtain the long distance pattern of 
a plane wave on the vortex, and then discuss its physical 
consequences. 

Wavefronts in the aeoustie ray approximation .— We 
now solve Eq. <11 using a ray, or WKB approximation. 


We assume the wave is propagating perpendicularly to 
the vortex axis. Eor this, we consider a solution of Eq. (§ 
independent of 2 : and of the form 

Uuj = and = ( 8 ) 

where Uk and Vk are (real) slowly varying amplitudes and 
S{r^0) the wave action. Putting this ansatz in Eq. 
and performing a gradient expansion we find that S is 
solution of the Hamilton-Jacobi equation 

{uj-v-VSf=F({VS)‘p. (9) 

In general this equation is solved by separation of vari¬ 
ables, i.e. assuming a single azimutal number. Since we 
want to consider an incoming plane wave, we shall instead 
obtain an approximate solution by treating the velocity 
field as a perturbation. Since ^ 0 at infinity, this will 
become valid far enough from the center of the vortex. 
As mentioned above, in the absence of velocity field, we 
simply have VS = k. In the presence of the vortex, 
we assume that the perturbation on the wave action is 
proportional to If 


VS = k -\- jff, ( 10 ) 

where 7 is position independent. Since ^ 0 at infinity, 
this equation means that there is an incoming plane wave 
of momentum k on the vortex. Assuming ylT is a small 
correction, we solve the Hamilton-Jacobi Eq. ^ pertur- 
batively. Expanding to first order, the derivative of F 
gives the Bogoliubov group velocity via Vg = kdkFjuj. 
Hence, we obtain the expression of 7 = — and deduce 
the wave action S{r,0) that governs the phase of WKB 
modes in Eq. ^ 

S{r, 0) = k ’ F — aO, ( 11 ) 


where 


a = 


kCw 


( 12 ) 


A key point is that the discontinuity in 0 must be the half 
line directed by k (we shall briefly return to that below). 
When the wave comes from the left, k - F = r cos(d) and 
d g] 0; 27r[. The key point of this result is that despite the 
fact that ^ 0 far from the vortex, there is a nonzero 
imprint on the profile of the modes, that is, they do not 
reduce to a plane wave. This is due to the fact that F falls 
off as 0 (l/r) at large r, something reminiscent to long 
range potentials in scattering theory m- The parameter 
a governs a long range effect induced by the velocity 
profile In addition, under the same approximations. 


^ The operator {iT-V} should be understood as (T • V + V • T‘)/2. 


^ Such a WKB method applied higher derivatives equations is ex¬ 
plain in great details in | 18| . see in particular Appendix A. 
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Figure 1. Hydrodynamic Aharonov-Bohm effect on a singly 
quan tize d vortex with w = +1: amplitude of SO or Sn of 
Eq. (16). Scale in the x and y directions is in units of the 
healing length. Red dots indicate the position of a vortex 
turning anti-clockwise, (a) = l/(47r) thus a ~ 1.6 and (b) 

= 1 thus a ~ 1. Depending on a of Eq. (12), one can 
observe a dislocation and / or a discontinuity in the forward 
direction of the plane wave. 


Eq. ^ gives a linear relation between the amplitudes 
Vk = DkUk, with 


Dk = (2^ky^lT^ - - l) . 


(13) 


Analogue Aharonov-Bohm effect in a BEC .— When 
using Eq. (11) for WKB modes ([^, we see that a plane 
wave scattered on a vortex is distorted by the second 
term of Eq. (11). This is the hydrodynamic analogue of 


the Aharonov-Bohm effect. When travelling around the 
vortex, for 0 from 0 to 27r, we accumulate a geometric 
phase of 271 a. In the original Aharonov-Bohm effect, the 
wave function is not a direct observable, and the elec¬ 
tromagnetic gauge symmetry only allows us to observe 
a modulo an integer [T]. Here, the perturbation field is 
a direct observable and extra features of the wave are 


described by Eq. (11). First, there are wave front dislo¬ 


cations governed by the integer part of a, i.e. there is [a] 
extra wave fronts coming out from the vortex core. This 
is very similar to the dislocation that one observes on the 
condensate wave function itself. Second, there is a line 
of discontinuity on ^ = 0 when a is not an integer. This 
feature cannot be seen on the condensate wave function, 
but only in the perturbations. 

Using the Bogoliubov dispersion relation, we compute 
the group velocity and use Eq. (§ to deduce the expres¬ 
sion of the Aharonov-Bohm parameter (p!^. 


a = 2^kw 


1 + 2C2fc2 • 


(14) 


It is instructive to explore this result in two opposite 
regimes: the long and short wavelength limit. 


a ‘ 


J 2w^k if ^k 1, 
\ w if ^k ^ 1. 


(15) 


In the long wavelength limit, the phase shift vanishes. 
This is because long wavelength excitations cannot be 


distinguished from the condensed part itself, i.e. they 
rotate the same way, and therefore, the relative phase 
disappears. For short wavelengths, the phase shift re¬ 
duces to 27r times the winding number of the vortex. 
This is due to the fact that wavelengths shorter than the 
healing length correspond to single particles, barely see¬ 
ing the condensate. Hence, they propagate without any 
extra phase, and the geometric phase of relative fluctu¬ 
ations is simply that of the condensate wave function. 
Additionally, this geometric phase affects both density 
and phase fluctuations of the condensate. Indeed, from 
Eq. 0, density fluctuations Sn are given by the real part 
of perturbations and phase fluctuations SO are given by 
the imaginary part. By choosing the phase convention 
such that Uk is real, we obtain 


SO 

Sn 

2no 


= (1 — Dk)Uk sin (^—ujt k ' 7 — aO^ , (16a) 
= (1 + Dk)Uk cos (^—ujt + k '7 — aO^ . (16b) 


Therefore, both observables show the same geometric 
phase shift, that is the analogue Aharonov-Bohm effect. 
The amplitudes of the perturbations are both governed 
by U/c. This relation between Sn and SO is a direct con¬ 
sequence of the continuity equation. 

Beyond WKB: the scattered wave. — Since 7 and its 
derivative vanish for r ^ oo, the WKB approximation 
becomes exact at long distances. In fact, WKB modes 
gives the incoming part of scattering states. A complete 
scattering state includes also a scattered wave and has 
the asymptotic form 


Akr 


(r, 0) ~ e*+A (0)^. (17) 

This is interpreted in a standard way m- The first term 
represents the incoming plane wave while the second one 
give the scattered wave whose amplitude fk{0) varies for 
different directions. Eq. 0 shows that the incoming 
part of stationary states does not reduce to a plane wave, 
but picks up the geometric phase governed by Eq. (11), 


something that was unnoticed in previous studies of scat¬ 
tering on hydrodynamic vortices [miinHii]. The second 
term of ( [l7| ) interfere with the incoming part to modify 
the wave profile. By comparing our Fig. with that of 
e.g. ID, we see that this mildly affects the main features 
described above. 

Moreover, when a is not an integer, the expression 0 
is only valid away from the line ^ = 0, where there is a 
discontinuity of the wave fronts (see Fig. (a)). When 
including the scattered part, the scattering amplitude di¬ 
verges in the forward direction, i.e. fk{0 -A 0) ^ oc. 
This shows that the WKB approximation breaks down in 
the forward direction. When considering exact solutions, 
the asymptotic (^ is valid only outside a small angu¬ 
lar region of width AO = 0{{kr)~^^‘^) where the mode 
smoothly interpolates between the WKB wave fronts on 
each side [T]. 
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Figure 2. Hydrodynamic Aharonov-Bohm effect on a hexag¬ 
onal vortex lattice: amplitude of SO or Sn on 25 vortices ar¬ 
ranged to a hexagonal lattice. Scales in the x and y direction 
is in units of the healing length, and vortices are separated by 
approximately 2^. Red vortices are turning anti-clockwise, 
and blue vortices clockwise. (a,b) = 1/5, (c,d) = 1, 

(e,f) — b. (Notice on the right panels the total circulation 

is one anti-clockwise vortex.) 


where Oj G]0,27r[ is the angle between 7 — r] and k. 
At large r, the Oj^s become the same, showing that the 
total geometric phase is given by the sum 27r aj . In 
Fig-i we plotted an example for a lattice consisting 
of 25 vortices with alined (left column) and alternating 
(right column) vortices. 

Outlook and summary .— The observation of solitonic 
vortex formation, as for example created by the Kibble- 
Zureck mechanism [23] , play an important role to further 
advance our understanding of quantum turbulence. The 
hydrodynamic Aharonov-Bohm effect presented above 
offers an alternative method to the homodyne detection 
method [iniiii] for quantum vorticity. Within the ho¬ 
modyne detection method the condensate is splitted into 
two parts. After a period of free expansion the two parts 
interfere and from the interference pattern one can di¬ 
rectly read off the dislocations /anti-dislocations of the 
fringes which correspond to the vortices in the initial 
condensate m- Based on the analysis above it is pos¬ 
sible to probe the vorticity, i.e. winding numbers and 
orientation of vortices, in BECs via the dislocation and 
discontinuity of phase- or density-fluctuations generated 
by a plane wave propagating through the condensate. 
The fact that the dislocation is reaching its maximum 
value — the number of new wave fronts emerging from 
a vortex is equal to its winding number — for ^ 1 
(i.e. high frequency excitations) implies that hydrody¬ 
namic Aharonov-Bohm can be utilized as a probe for 
quantum vorticity. The experimental advantages, when 
compared to the homodyne technique, are that the hy¬ 
drodynamic Aharonov-Bohm effect can be observed di¬ 
rectly as matter waves in the condensate. 

The detection of the hydrodynamic Aharonov-Bohm 
effect seems to be within experimental reach. Recent 
works have shown very promising results, in both the 
capacity to excite [24| and image [25l [26] plane waves in 
trapped condensates. When probing quantum vorticity 
with sub-healing wavelengths, deviations from the 
predicted pattern are expected, that might reveal details 
of the core of the quantum vortex. 


Multiple vortices .— Although our results are valid for 
any vortex, in a condensate, only single quantized vor¬ 
tices (I re I = 1) are stable [12]. When the condensate is 
spinning faster, higher vortices decay into several single 
quantized vortices. One can however obtain a stable con¬ 
densate with a lattice of vortices. The linear character 


of the long range effect described by Eq. (10) shows that 


the hydrodynamic Aharonov-Bohm effect cumulates and 
is only sensitive to the total circulation. More precisely, 
we consider a lattice of vortices located at fj and of asso¬ 
ciated Aharonov-Bohm parameter aj. The phase of the 
Aharonov-Bohm mode ^ reads 


S{r,0) = I -r 


(18) 
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